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o : 

CSl i Abstract. A notion of monitored recurrence for discrete-time quantum processes 

' was recently introduced in |15| taking the initial state as an absorbing one. We 

, extend this notion of monitored recurrence to absorbing subspaccs of arbitrary finite 

pLn ' dimension. 

, The generating function approach leads to a connection with the well-known 

' theory of operator-valued Schur functions. This is the cornerstone of a spectral 

characterization of subspace recurrence that generalizes some of the main results in 
(-H ' |15j . The spectral decomposition of the unitary step operator driving the evolution 

^ |. yields a spectral measure, which we project onto the subspace to obtain a new 

I ' . ' spectral measure that is purely singular iff the subspace is recurrent, and consists 

Ch . of a pure point spectrum with a finite number of masses precisely when all states 

^ ' in the subspace have a finite expected return time. 

This notion of subspace recurrence also links the concept of expected return time 
to an Aharonov-Anandan phase that, in contrast to the case of state recurrence, can 
be non-integer. Even more surprising is the fact that averaging such geometrical 
" ^ ' phases over the absorbing subspace yields an integer with a topological meaning, 

, so that the averaged expected return time is always a rational number. Moreover, 

QQ ' state recurrence can occasionally give higher return probabilities than subspace 

, recurrence, a fact that reveals once more the counterintuitive behavior of quantum 

systems. 

^SJ , All these phenomena are illustrated with explicit examples, including as a natural 

' application the analysis of site recurrence for coined walks. 

m ■ 

• • ' 1. Introduction 

^> : 

^ ■ The study of return properties of discrete-time random processes goes back at 

^ . least to G. Polya in 1921 [21], who proved that only in dimension one and two is 

the simplest unbiased random walk recurrent, returning to the starting point with 
probability one. Since then, the notion of recurrence has been key in the theory of 
classical random processes. 

Concerning the analog of recurrence for discrete-time quantum systems, two dif- 
ferent approaches have been recently proposed. In 2008, M. Stefahak, I. Jex and 
T. Kiss [211 E21 122] resort to an ensemble of identically prepared independent systems 
that are discarded after the measurement to avoid changing the dynamics due to the 
collapse. A more recent proposal, from A.H. Werner, R.F. Werner and two of the 
authors [15], treats the perturbation of the evolution due to the monitoring as an 
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essential ingredient in the definition of quantum recurrence. Tliis second definition 
of recurrence is based on performing a projective measurement after each step of the 
quantum evolution, an idea that has some precursors in the context of quantum walks 
with absorbing boundary conditions [H |3] . 

The ideas developed in [15] deal with the probability of returning to the initial 
state, which corresponds to having an absorbing one-dimensional subspace generated 
by this initial state. A natural generalization, which is the aim of this paper, is 
to consider an absorbing subspace of arbitrary finite dimension which includes the 
initial state. This is the case, for instance, when we are interested in the return of 
a coined walk to the initial site, regardless of the internal state. This generalization 
also allows us to ask about the probability of returning to a finite set of sites from 
which an initial state was chosen. We will use the terminology "state recurrence" to 
refer to the notion of recurrence defined in [15], and "subspace recurrence" for the 
extension of this concept developed in the present paper. 

The spectral characterization of subspace recurrence is the principal objective of 
this paper, which emerges from the identification of a key object that best codifies 
subspace recurrence, namely the operator-valued Schur function of the spectral mea- 
sure of the absorbing subspace. This connects the subject of quantum recurrence with 
well-known areas of complex analysis: the theory of matrix Schur functions [TU |5] 
and the theory of orthogonal polynomials on the unit circle [26| [27]. See also [13] for 
a more detailed discussion of the matrix version of these orthogonal polynomials and 
its relation with matrix Schur functions. 

It turns out that the recurrent subspaces - those whose states all return with 
probability one to the subspace - are characterized by inner Schur functions. This 
identifies the recurrent subspaces as those contained in the singular subspace of the 
unitary step operator driving the evolution. Furthermore, the recurrent subspaces 
whose states all possess a finite expected return time are those whose Schur functions 
are rational inner. Such subspaces are those contained in finite sums of eigenspaces 
of the unitary step. These results are the natural extensions of Theorems 1 and 2 in 

m- 

We stress that the aim of the paper is not limited to the generalization to subspace 
recurrence of the main results already known for state recurrence. The paper also 
shows important differences and unexpected relations between subspace and state 
recurrence. For instance, in contrast to state recurrence, the existence of states that 
return with probability one to a subspace of finite dimension greater than one does 
not require the existence of a singular subspace for the unitary step. This is true even 
if the state returns to the subspace in a finite expected time. 

Besides, in the case of state recurrence, the expected return time is always infinite 
or an integer. This is because the expected return time has a topological interpre- 
tation as the degree of the related inner Schur function |T5]. These results are no 
longer true for subspace recurrence because the expected return time of a state is 
related to a geometrical (rather than a topological) invariant in this case, namely the 
Aharonov-Anandan geometrical phase of the curve obtained by applying the (oper- 
ator) boundary values of the Schur function on the unit circle to the state. As a 
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consequence, the expected return time to a subspace can be any real number not 
smaller than one. 

However, a remarkable thing occurs when averaging the expected return time over 
the absorbing subspace: such an average is always infinite or a rational number. 
Actually, the appearance of integer quantities in the finite case has again a topological 
explanation: the numerator of the rational average is the degree of the determinant 
of the corresponding Schur function, while the denominator is simply the dimension 
of the subspace. This integer numerator may also be expressed in terms of invariant 
subspaces for the unitary step (as the dimension of the minimal invariant subspace 
containing the absorbing subspace), and in terms of the spectral measure of the 
subspace (as the sum of the ranks of all the mass points). 

We should also highlight a quantum paradox that arises from the comparison be- 
tween state and subspace recurrence: the probability of returning to a state can be 
bigger than to a subspace containing such a state. Indeed, a state can return to itself 
in a finite expected time without returning with probability one to the subspace (see 
Example 12. 3p . Thus, in contrast to the classical expectation, the return probability 
can decrease when the dimension of the absorbing subspace increases. Examples of 
this for one- and two-dimensional coined quantum walks will be given in Section H] 
(see Figs.EHU). 

The use of Schur function methods to the study of quantum recurrence is not limited 
to the spectral characterization of recurrence. For instance, the return probability to 
a subspace, as well as the expected return time, can be computed using the Schur 
function of the corresponding spectral measure (see Section H]). We find that many 
techniques and results borrowed from the theory of orthogonal polynomials on the 
unit circle are invaluable tools for computing this Schur function (see Appendix [B|) . 

Regarding the structure of the paper. Section [2] develops the notion of subspace 
recurrence, including its spectral characterization. The expected time for the return 
of a state to a subspace is discussed in Section [3l which ends with a spectral charac- 
terization of the subspaces whose states all return in a finite expected time. Section |4] 
deals with applications of the previous results to the study of site recurrence for coined 
quantum walks. Finally, Appendix |X] includes the proofs of some technical results of 
interest for Section [3l while Appendix |B] gives the details of the calculation of the 
matrix Schur function for a site in a one-dimensional coined walk, which is used in 
the applications of Section HI 

2. Subspace recurrence 

We will consider a quantum discrete-time evolution driven by a unitary step op- 
erator f/ on a complex separable Hilbert state space "H with inner product (-I-) and 
norm || ■ ||. The rays of "H, which are represented up to phase factors by vectors on 
the unit sphere S-u = {'ip ^ 'H : \\ip\\ =1}, are the quantum states that evolve after n 
steps via ip t/""?/'. 

Given a subspace V of Ti and a state ip & Sy = {'ip & V : \\ip\\ = 1}, we can ask 
about the probability of returning to the subspace V if the evolution starts at the 
state ip- Following the ideas developed in [13], a monitored procedure to define such 
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a probability consists in introducing measurements to check the return of the system 
after every step. 

More precisely, starting at ip = ipo, any step ipn Uipn is followed by a measure- 
ment of the orthogonal projection P = Py onto V . When this projective measurement 
succeeds, the experiment ends. Otherwise the system is projected to V-^ so that it is 
left in the state V'n+i £ Sq-i obtained by normalizing (/ — P)Uipn- 

If the system was not found in the subspace V during the first n steps, after the 
n-th step it will be in a state proportional to U^if) with U = {I — P)U. According to 
the standard rules of Quantum Mechanics, the probability of this survival is = 
||f/"?/'|p, and we will call it the n-step V -survival probability of ip. 

It follows that the probability that the system returns to the subspace V for the 
first time in the n-th step comes from multiplying the n — 1-step survival probability 
||fjn-i^||2 i^-y ^Yie probability UPt/^/^n-iH^ to find the system in the subspace V con- 
ditioned on such a survival. Therefore, the probability we are searching for is given 
by WPULf^'^ipW"^, and we will refer to it as the n-step first V -return probability of ijj. 

We emphasize the difference between the first V^-return probabilities and the quan- 
tities llPt/"?/'!!^, which also represent the probability of finding the system in the 
subspace V, but performing the projective measurement only in the n-th step, i.e. 
without any intermediate monitoring of the process. To distinguish both kind of 
probabilities we will refer to HPf/"?/'!!^ as the n-step V -return probability of tp. 

First V^-return probabilities and V^-return probabilities can be written as Ha.n^/'lp 
and ll/Li^V'p respectively, where a„ and /x„ are the operators on V given by 



a„ = PUU^'' P n-step first V -return amplitude operator, 
/x„ = PU^P n-step V -return amplitude operator. 



We will use boldface characters to denote those operators that, like a„ and should 
be considered as operators on V . Note that a„ and /x^, apart from ra, depend (only) 
on the unitary step U and the subspace V . If it is convenient to indicate explicitly this 
dependence for any of the mathematical objects that arise, we will use a superscript 
and a subscript, writing for instance a^y and ^i^y- 

The matrix elements of the above operators also provide transition amplitudes of 
interest. For any states ip,(t> E Sy, the squared modulus of the amplitude {(f)\ani^) 
gives the probability of the n-step transition ip (p without hitting V at intermediate 
steps, while the squared modulus of the amplitude {(plfi^^ip) yields the corresponding 
n-step transition probability without intermediate monitoring. 

To avoid overcounting probabilities, the total probability of a transition - the prob- 
ability that it will take place eventually, without any reference to when it takes place 
- should be defined as the sum over the number of steps n of the n-step first time 
transition probabilities. This leads to the following total probabilities for any states 
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e Sv, 

R{ip) = ll'^'n^ll^ V -return probability of ip, 

(2) 

Prob(?/', 0) = \{(l)\anip)\^ ip ^ (j) probability when returning to V . 

n>l 

The second quantity measures the relative hkehhood of landing on one state in V 
compared to another. Note that R{il)) = '^i.PToh{ip,(f)k) for any orthonormal basis 
(pk of V. 

Quantum physicists, convinced by their expertise on the consistency of Born's 
probabilistic interpretation of Quantum Mechanics, would say that the probabilistic 
nature of the above quantities ensures that J2n>i ll'^n'^lP J2n>i \{4>\f^n'ip)\'^ are 
not greater than one. Nevertheless, this can be easily proved, showing that the 
probabilistic interpretation of the above series does not lead to a contradiction. 

The orthogonal decomposition UU"'~^P = + tf^P yields 

||a„V>||2 = ||f/"-Vf - ^"V'f, i^eV, (3) 
thus the V^-return probability of a state ^/^ G is given by 

R{il)) = 1 - lim . 

This identity not only proves that R{ip) < 1, but also that Sn{ip) = \\U""ip\\'^ converges 
to a number s{ip) = ^^^n-^oo llf^""^!!^ < 1 that provides the total V -survival probability 
of ip, i.e. the probability that the state ip is never recaptured by the subspace V. 
The convergence of the n-step survival probabilities Sni'ip) to a number not greater 
than one can also be proved directly by checking that ||f/"?/'|| is non-increasing and 
bounded by one. Finally, the fact that |(0|a„?/;)| < ||a,.„^/^|| for any (p E Sy yields 
Prob(V',0) < 1. 

Example 2.1. Let H = ©£^(Z) with an orthonormal basis whose states we denote 
by It)) \i) and \x), X G Z. We define a unitary operator U on H as the shift 
U\x) = |x + 1) on the subspace spanned by \x), x G Z, together with f/|t) = ||), 
U\i) = It)- 

Consider the subspace V = spa.ia{tp, </>}, where 

^ = n), = -^(it) + io)). 

Then, [7^ = i(||) - |0)) and t/"V^ = -i|n - 1) for n > 2, while = ^|n) for 
n > 1. Therefore, 

r.i?/) = -^A) I n„ih — li/i 

a„, = for n > 3. 

For any ipa,i3 = Oiifj + f3(f) E Sy we find that R{ipa,i3) = | + \ which reaches its 
maximum 3/4 and minimum 1/2 at -j/^ and 0, respectively. The transition probabilities 
when returning to V between these two extreme states are Prob(^/', 0) = Prob(0, = 
1/2. ■ 
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2.1. Generating functions. The inequahties ||At,J|, ||an|| < 1 ensure that the series 

f^(^) = Yl ^^^^ = Yl 1^1 < 1' 

n>0 n>l 

are absolutely convergent and hence define operator-valued functions on the unit disk. 
We will refer to jl{z) and a{z) as the l^-return and first ^-return generating functions, 
respectively. 

These generating functions are related to the Green functions of the operators U 
and U, 

G{z) = {I - zU)-\ G{z) = {I-zUy\ 
through the relations 

^i{z) = PG{z)P, a{z) = zPUG{z)P. 

The above relations, together with the resolvent identities 

G{z) - G{z) = zG{z)PUG{z) = zG{z)PUG{z), 

and the equality PG{z)P = P, which follows from G{z) — I = zUG{z), lead to 

fi{z) -P = fi{z)a{z) = zPUG{z)P. 

Expanding G{z)P = G{z){P + zPUG{z)P) = G{z)Pfi{z) in the last term yields the 
renewal equations 

fi{z)a{z) = a{z)fi{z) = fi{z) - P, 

which are the extension to dim > 1 of the renewal equations already found in [15] 
for diml^ = 1. In terms of the Taylor coefficients, these equations become 

n—l 1 

H^ = Un + Y t^k^n~k = an + Y "-n-kfJ-k- 

k=l k=l 

In the classical case, an identical formula expresses the probability of returning in 
n steps in terms of all the ways it can return for the first time in n or fewer steps. 
However, in the quantum case, a„ and fi^ are (operator) amplitudes rather than 
probabilities. Nevertheless, both renewal equations hold despite the fact that the 
operator coefficients /z„ and a„ do not commute in general. In any case, they allow 
us to compute the first ^-return amplitudes a„ from the l^-return amphtudes 
which are usually more accesible. 

Since we consider the coefficients a„ as operators on V, the generating functions 
fi{z), a{z) become functions with values in operators on V. Thus the term P in the 
renewal equations must be understood as the identity ly in V, leading to 

{Iv - a{z))fi{z) = ly. 

This shows that fi{z) and ly — o,{z) are invertible on the unit disk with 

= i^v - a{z))-^, a{z) = Iv- K^Y^ ■ 
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These relations permit us to express the l^-return amphtudes purely in terms of the 
first ^-return amplitudes and vice-versa, 



ni+n2=n ni+n2+n3=n 



an = Hn- f^n^f^n2 + /^^iMnzMna + i'^T Vl- 

ni,n2>l ni,n2,n3>l 
ni-j-n2=n ni+n2+n3=n 

We can arrive at a spectral interpretation of the above generating functions from 
the spectral decomposition of the unitary operator U, 

U = [ \E{d\), I E{d\) = I, 

which expresses U as an integral over the unit circle = {z ^ C : \z\ = 1} with 
respect to the projector- valued spectral measure E{dX). This allows us to rewrite its 
Green function as 

E{d\) 
l-\z' 



51 



G{z) 

The relation between G{z) and Ij,{z) gives 

A(^) = / f^idX) = PEid\)P, (4) 

Jgi i — Az 

where fi{d\), which must be considered as a measure with values in the space of 
operators on V, will be called the spectral measure of the subspace V. The measure 
l-i{d\) is supported on a subset of the unit circle and, since E{dX) is a resolution of 
the identity, J^i i-i{dX) = Jy. For convenience, we will omit the domain of integration 
with respect to fi{dX) because it will always be S^. 

The equality (jlj) identifies fj,{z) as the operator-valued Stieltjes function of the 
measure ^{dX). When V is of finite dimension, identifying an operator on V with a 
finite matrix allows us to use the results of the theory of matrix-valued Caratheodory 
and Schur functions [HIE]. In what follows we will assume that dimy < 00. 

With this assumption, the operator-valued Caratheodory and Schur functions of 
the measure ^j,{dX) are the analytic functions on the unit disk defined by 

F{z) = t^idX), f{z) = -{F{z) - Iv){F{z) + Iy)-\ \z\ < 1. (5) 

J X - z z 

We will also refer to F{z) and f{z) as the Caratheodory and Schur functions of the 
subspace V. Evidently, F(0) = ly, ReF{z) := + F{z)^) > and F^{z) = 

2fi(z) — Iv, where F\z) = F(zy is the analytic function on the unit disk whose 
Taylor coefficients are the adjoints of those of F{z). 

The relation between Caratheodory and Schur functions gives 

{Iv-zf{z)){F{z) + Iv) = 2Iv, 

ReF{z) = (ly - Zf{z)^)-\ly - \z\'fizyfiz)){Iy - zf{z))-\ (6) 

ly - \z\'f{z)^f{z) = 4(F(^)t + ly)-^ Re F{z){F{z) + ly)-' > 0, 
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which shows that f[zYf{z) < ly due to Schwarz's lemma. This inequahty is equiv- 
alent to f{z)f{zy < ly, and can be also written as < 1. If f{z) is not 
constant, the maximum principle implies that these inequalities are strict inside the 
unit disk. Constant Schur functions appear for instance when V is ^/-invariant (see 
below). 

From (|6]), we obtain fi{z) = (ly — 2/^(2;))"^. Thus, the renewal equation can be 
rewritten as 

a{z) = zf{z). 

This means that the first V -return generating function is the Schur function ofV, up 
to conjugation and multiplication by z. In other words, the first l^-return amplitudes 
a„ are the adjoints of the Taylor coefficients for the Schur function f{z) of the sub- 
space V. In the case dim^ = 1, this relation between the first ^-return generating 
function and the Schur function of V was reported for the first time in |15j . 

When V is [/-invariant, P commutes with [/ in ([T]) and a{z) = zPUP. This 
identity refiects the obvious fact that any state of V returns to V in one step, and 
ai = U \ V is the restriction of U to V. We also see that f{z) = PWP is constant 
in this case. 

The boundary behavior of Flz) as \z\ — t- 1 provides relevant information about the 
measure: if fi{d\) = w{9)^ + fj,^^^^{dX), A = e*^, is its Lebesgue decomposition, then 
w{9) = limj,_>.i- Re-F(re*^) for a.e. 9; the singular part fi^^^^{d\) is concentrated on 
the points z E such that lim^^i- Tr Re -F(r2;) = 00; and the mass points z E 
are characterized by a non-zero value of fJ^dz}) = hm^^i- ^^F{rz). These results 
follow from the analogous ones for the scalar case using the polarization identity for 
an operator T on a complex Hilbert space, {(j)\Tip) = ^ Ylk=o ik{4> + + ^V^)) ; 

see [T3] . 

Like Caratheodory functions, Schur functions can be defined a.e. on the unit circle 
by their radial limits, /(e*^) = lim^.^!- /(re*^). The following proposition gives a 
central result for the analysis of quantum recurrence. 

Proposition 2.2. The measure fi{dX) is purely singular exactly when its Schur func- 
tion f{z) is inner, i.e. when /(e'^) is unitary for a.e. 9. In particular, ^{dX) is 
finitely supported - it is a sum of finitely many mass points - if and only if f{z) is a 
rational inner function. 

Proof. Taking boundary values in ([6]) we find that Rei^(e*^) = is equivalent to 
= ly- Thus, up to a set of zero Lebesgue measure, w{9) = if and 
only if /(e*^) is unitary. When ^{dX) is finitely supported, F{z) and f{z) are obvi- 
ously rational. Conversely, a rational Schur function f{z) corresponds to a rational 
Caratheodory function F{z). Since F{z) is rational, /Ligjj^g((iA) must be concentrated 
on the set of finitely many poles of F{z) that lie on the unit circle. Since f{z) is 
inner, w{9) = a.e. and fi{dX) = /-i^i^^^dX) is finitely supported. □ 

2.2. V-return probability and V-recurrence. The previous result is the key for 

the use of techniques of complex analysis in the study of quantum recurrence. A 
first consequence is that the generating function a{z) is also a Schur function with 
radial boundary values a(e*^) = e*^/(e~*^)^ a.e. on the unit circle. The transition 
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probabilities related to the subspace V can be expressed in terms of these boundary 
values as 

where R{ip) and Prob(?/',0) were defined in above. Therefore, the V -return prob- 
ability R{ip) = {ip\Rip) is a quadratic form in V given by the operator R, which can 
be understood as the squared 'right operator norm' of the Schur function f{z). 

We will refer to R as the V -return probability operator, which obviously satisfies 
< -R < /y in agreement with the fact that {ip\Rip) must be a probability for any 
ip G Sv- As a consequence, when ip runs over the unit sphere Sy, the possible values 
of R{ip) cover the entire interval between the minimum and maximum eigenvalues 
of R, and these eigenvalues are the minimum and maximum l^-return probabilities. 
If these eigenvalues are distinct, i.e. R{ip) is not constant on Sy, the corresponding 
eigenspaces are orthogonal to each other and the states of V with maximum and 
minimum ^-return probabilities are the unit spheres of these two eigenspaces. 

The mean value of the ^-return probability is the unbiased average of R{ip) over 
the states of V, namely 

/ RWd^, 

JSv 

with dip the constant probability measure on Sy. Averaging the diagonal coefficients 
(iplTip) of any normal operator T on over Sy is equivalent to averaging over any 
orthonormal basis of V. Indeed, if (pk is an orthonormal eigenvector basis of T and 
Afc are the corresponding eigenvalues, 

£WT,P) * = X: A. £ mWd,p = = J^^. (7) 

where we have used that 

Applying this to the positive semidefinite operator R we find that 

where \\A\\f = TtIA'^A) is the Frobenius norm. This result states that the mean 
V -return probability equals the arithmetic mean value of the eigenvalues of the V- 
return probability operator R, counting multiplicity. Equivalently, the mean l^-return 
probability comes from averaging over the unit circle the arithmetic mean of the 
squared singular values of the Schur function f{z) of V , counting multiplicity. 

We say that a state ip & Sy is V -recurrent if it returns to V with probability one, i.e. 
R^ip) = 1. In particular, the recurrence of a state ip, as defined in corresponds 
in our context to the special case that V is the subspace spanned by ip. One of the 
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main results in [15] is that the set of recurrent states is the unit sphere of the singular 
subspace of U. As we will see below, when dim V > 1, the l^-recurrent states continue 
to possess a subspace structure (up to normalization), but it need not be the singular 
subspace of U. 

Due to the identity R{ip) = {4-'\Rtp) and the inequalities < -R < Jy, the existence 
of ^-recurrent states is equivalent to the condition = 1. This means that 1 is 
an eigenvalue of R and the unit sphere of the corresponding eigenspace is the set of 
V^-recurrent states. Similar arguments show that a state ip is ^-recurrent if and only 
if ||a(e*^)^/;|| = 1 for a.e. 9, which means that f{e^^)f{e'^^yip = ip except on a set of 
zero Lebesgue measure. Therefore, the existence of l^-recurrent states implies that 1 
is a singular value of /(e*^) a.e. The converse is not true because the corresponding 
eigenvectors of /(e*^)/(e*^)^ will generally depend on 6. 

If fj,{dX) = w{9)^ + /^gi^g((iA), we know that the weight is given for a.e. A = e*^ by 
w{6) = Rei^(A). The existence of ^-recurrent states implies that ||/(2)|| = 1 a.e. in 
S^, which, in view of ([6]), is equivalent to saying that the operator w{6) is singular 
for a.e. 6. Hence, the presence of ^-recurrent states requires an a.e. singular weight 
for the absolutely continuous part of the spectral measure of V. The above comments 
suggest that the converse is not true. (This is confirmed in Example 12.31 below). 

Thinking classically, we expect the l^-return probability of a state to be non- 
decreasing when enlarging the subspace V. However, in quantum recurrence, mod- 
ifying the target subspace V also changes the monitoring of the process. This has 
the effect of dramatically altering not only the possible return paths to V, but also 
the way in which each path contributes to the total l^-return probability. The reason 
for this is that Quantum Mechanics adds complex amplitudes or positive probabili- 
ties depending on whether the paths have the same target state or orthogonal target 
states. This results in the absence of monotonicity with respect to the subspace V 
for the l^-return probability. In particular, the recurrence of a state ip is in general 
not easily connected to its V^-recurrence for a subspace V that contains ip. 

We will see that the above effects already appear in simple models like the one 
introduced in Example 12.11 We will also use this example to illustrate the generating 
function approach to quantum recurrence. Although this approach could be avoided 
for such a toy model, the analysis below will clarify in a simple setting some aspects 
that will reappear in the applications discussed in Section HI where Schur functions 
become an indispensable tool for the analysis of recurrence. 

Example 2.3. Let us compute the generating functions for the subspace V = spa.n{ip, 0} 
of Example 12.11 Identifying any operator on V with its matrix representation in the 
basis {ip,(p}, and denoting by I2 the 2x2 identity matrix, we have the l^-return 
amplitudes 
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which lead to the generating functions 
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This gives the l^-return probabihty matrix 

-.27r 

R 





I 

i 



in agreement with the V^-return probabihty obtained in Example 12. 



The corresponding Caratheodory and Schur functions 





1 


1, 




= 2' 


4 



Fiz) 



■ ' l-z- 

yield the spectral measure of V 

de 1 
2^ ^ 2 



1 + ^2 



n{dX) 




i 



1 



V2z 
1 



1 

2 



£ J_ 

2 V2 

4^ 



V2 



5(A-1) + - 



-1 



-1 

2 



<5(A 



where A = e . Note that the weight w{6) is singular despite the absence of V- 
recurrent states. This occurs because 1 is a singular value of /(e*^) but the corre- 
sponding eigenvector of /(e*^)/(e*^)^ is 6'-dependent. 

We can compare state and subspace recurrence by computing the return probability 
Ra^l3 of the state 4'a,i3- The scalar generating functions 



P'aAz) = ( 

da,fi{z) = 1 

allow us to calculate 



a I3)(i{z) 
1 



1 - 



j2 + ^2 Re (a/3); 



(1 



2) z2 + y2Re(a/3); 



2tt 



j2 + y2Re(a/3); 



Ra,i 



JS\\2 



de_ _ i-\w 

1/3 1 2 (see Figure [I]). We find that Rc,,^ > 



which should be compared to R{ipa,p) = 4 : 
R{'^a,/3) when |/3p < |(5 — vTf) ~ 0.44, which shows that the return probability to a 
state can be bigger than the return probability to a subspace containing such a state. 
In fact, ip & V is a recurrent state while V has no V -recurrent states! Furthermore, 
ijj is not l^-recurrent despite the fact that the expected time for the return of ijj to 
itself is finite and equal to 2 (see Section |3]). ■ 

The previous example shows explicitly that a subspace V with recurrent states 
need not have any ^-recurrent states. On the other hand, it is easy to check that 
the subspace V = span{|0), |1)} in Example 12.11 has no recurrent states while |0) is 
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Figure 1. State versus subspace recurrence in Example 12. 3[ 

The return probability for tpa^/s = aip + P(f) ^ function of \(3\. The 
return probability -Rq,^ to the state (dashed blue line) is higher than 
the return probability R{ipa,i3) to the subspace V = spanj?/;, 0} (solid 

red line) when \(3\ < J \{h - y/Vl) ^ 0.66. 



recurrent. Thus, a subspace with ^-recurrent states need not have any recurrent 
states. Therefore, state recurrence and subspace recurrence seem to be unrelated 
notions. Nevertheless, we will see that, surprisingly, demanding all the states of V to 
be recurrent or V^-recurrent are equivalent requirements. 

Theorem 2.4. Let U be a unitary step and V a finite- dimensional subspace with 
spectral measure fi{dX) and Schur function f{z). Then, the following statements are 
equivalent: 

(1) All the states ofV are V -recurrent. 

(2) All the states of V are recurrent. 

(3) fi{dX) is purely singular. 

(4) f{z) is inner. 

Proof. Bearing in mind that /(e*^)/(e*^)^ < Jy for a.e. 6 we find that 
R{tlj) = eV <^ R = ly f{z) is inner <^ 

■x^ fj,{dX) is purely singular. 

On the other hand, the spectral measure of a state is given by fi^{d\) = {ip\E{d\)ip) , 
where E{dX) is the spectral measure of U. If iJj ^ Sy, then clearly fi^{dX) = 
{ip\fi{dX)ip) . Hence, the polarization identity shows that fj,{dX) is purely singular 
if and only if ^^{dX) is purely singular for any ip G Sy- This is equivalent to stating 
that all the states of V are recurrent. □ 
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We will say that a subspace V is recurrent if every state if) ^ Sy returns to V 
with probability one, i.e. if every state ip E Sy is ^-recurrent. This notion obviously 
coincides with the mere existence of ^-recurrent states for the case dim V = 1 defining 
state recurrence; otherwise, the presence of l^-recurrent states does not necessarily 
mean that they exhaust all of Sy- 

The above theorem is the generalization of [T5| Theorem 1] to subspace recurrence. 
It states that the finite-dimensional recurrent subspaces are exactly those whose states 
are all recurrent. As a consequence, any subspace of a finite-dimensional recurrent 
subspace is recurrent too. Another way to restate the previous result is to say that the 
finite-dimensional recurrent subspaces are those contained in the singular subspace 
of the unitary step U. 

We stress that despite the close analogy between Theorem 12 . 4 1 and p!5| Theorem 1], 
new situations appear when extending the concept of state recurrence to subspaces, 
even if one restricts oneself to a finite-dimensional subspace V. In contrast to the case 
of state recurrence, the presence of l^-recurrent states does not imply the existence 
of a singular subspace for U since such l^-recurrent states need not cover all of Sy- 
This is clearly shown, for instance, by the shift U\x) = |x + 1), x G Z, in since 
the state |0) is ^-recurrent for V = span{|0), |1)}, but U has no singular subspace. 
Moreover, the existence of a singular part for the measure ^{dX) of a subspace V 
does not ensure the presence of ^-recurrent states. Example 12.31 illustrates this fact 
because V = spanj?/', 0} is free of l^-recurrent states, while the related measure i-i{dX) 
has two mass points, but is not purely singular due to the absolutely continuous part. 
Additional examples of these phenomena will appear in Section HJ 

3. Expected return time to a subspace 
The expected l^-return time of a state ip G Sy, defined by means of 

n>l 

makes sense whenever ip is l^-recurrent, i.e. Yln>i ll'^n'^P ~ Otherwise, T{ip) is 
naturally taken to be infinite. 

Like the l^- return probability R{ip), the expected l^- return time of any state if) eV 
can be expressed in terms of the ^-survival probabilities = llf^^V^P- Using once 

again the relation (|3]) we find that 

N N-1 
n=l n=0 

Therefore, the series on the right hand side must be divergent whenever T{ip) = 
oo. On the other hand. Lemma IA.3I of Appendix [B] shows that r(V') < oo implies 
lim^^oo NWU^ipW"^ = 0. Hence, the equality 

n>0 
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always holds, which means that the expected V -return time of a V -recurrent state is 
equal to the sum of its n-step V -survival probabilities. 

We can also express T{ip) in terms of the Schur function f{z) of V. For any V- 
recurrent state ip we can write 



In particular, if V is recurrent and its inner Schur function f{z) has an analytic 
extension to a neighborhood of the closed unit disk, e.g. if f{z) is a rational inner 
function, then we can write 



where ^l>{0), 9 G [0, 2tt], traces out a closed curve on the sphere Sy due to the unitarity 
of d(e*^). This simple result has a nice interpretation since it relates T{ijj) to a kind of 
Berry's geometrical phase [11 [25]. More precisely, the expected V -return time of 
a state ip & Sy is —1/2-11 times the Aharonov-Anandan phase [1] associated with 
the loop a{e'^)ijj: Sy 

In the case of state recurrence, [I5] proves that the states ip with a finite expected 
return time are characterized by a finitely supported spectral measure fi^,{d\), thus 
by a rational inner Schur function f^{z). Further, [15] also finds that r (■?/') must be 
a positive integer whenever it is finite because of its topological meaning: t{iI)) is 
the winding number of a^(e*^): — )■ S^, where d^{z) = zf^{z) is the first return 
generating function of if). 

In contrast to a winding number, the Aharonov-Anandan phase is not necessarily 
an integer because it reflects a geometric rather than a topological property of a closed 
curve. The expression ([H]) for t{iI)) is reparametrization invariant, and changes by an 
integer under closed gauge transformations ip{9) — >■ ip{9) = e*'^^^V(^)) fpi^ir) = 
ip{0). This means that t{i/j) is a geometric property of the unparametrized image of 
ip{9) in Sv, while e*^'^'^'^'^) is a geometric property of the corresponding closed curve in 
the projective space of rays of Sy whose elements are the true physical states of V. 
In fancier language, Sy is a fiber bundle over such a projective space with structure 
group S^, and e~^'^'^'^^'^^ is the holonomy transformation associated with the usual 
connection given by the parallel transport defined by {ip{t)\dtip(t)) = [1]. 

As a consequence, we cannot expect for r(?/^) to be an integer for subspaces V 
of dimension greater than one. The following example shows explicitly that the 
expected \^-return time can be non-integer valued, a fact that will also be illustrated 
by the applications discussed in Section HI Like Example 12. 3[ the following example 
demonstrates the generating function approach to computing the expected l^-return 
time on a simple problem where other methods are available. 

Example 3.1. Given a three-dimensional Hilbert space "H with orthonormal basis 



{|0), |1), |2)}, let U be the cyclic shift U\0) = |1), U\l) = |2), U\2) = |0). Any 





d9 
2^' 
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subspace of H is recurrent in this case because the corresponding spectral measure 
is necessarily finitely supported. Therefore, the expected ^-return time makes sense 
for any state V'o,/? = tt|0) + (3\l) G Sy of the subspace V = span{|0), |1)}. Then, 
Uipa 13 = /9|2) and t/'^Va /? = for n > 2, while the first V^-return amplitudes are given 
by 

ai^a,/3 = a|l), o,2ipa,i3 = /3\0), a„ = for n>3. 
We conclude that 

which is an integer only for the states |0) and |1). 

We can recover this result using the generating function approach. Identifying any 
operator on V with its matrix representation in the basis {|0), |1)}, the l^-return 
amplitudes 













/^3n+l — \l 


or 






n 



t^zn — fJ'3n+i — 1 1 n ) ' ^3n+2 " ( n n ) ' n>0 
yield the generating functions 



Therefore, we can recover the previous expected V^-return time from 

r(^„,,) = (a ^) r , r = a(e^^)ta,a(e^^)|^ = °) . 

We can also compute the average of the expected l^-return time over the sphere 
of states Sv- If <j{dad/3) is the constant probability measure on the 3-sphere = 

{{a,P)eC':\a\' + \P\' = l}, 



T{^)dij= [ il + \P\')aidadP) = l, 



'Sv JS^ 

which is not an integer but a rational number. ■ 

The above example suggests that, although the geometrical phases giving the ex- 
pected V^-return time can be arbitrary, their average over the states of V is always a 
rational number. We will see that this is indeed the case, but before doing that we 
will express such an average in a more convenient way. 

Using the monotone convergence theorem for dtp and the relation we find 
that the averaged expected l^-return time is 



/ T{ip)dip = -f nWanipfdip = / y^n||a„||| 

^— lim f \i(a{re'ydea{re'^)) — , 
im V r--s>i- Jq 2'iri 

l + -i-lim r TT{def{re^')f{rey) 
dim V r^i- Jq 



dim 

1 .....x d9_ 

where we assume that the subspace V is recurrent. 
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Apart from its own interest, this average is a useful tool to characterize those 
recurrent subspaces V with T{ip) < oo for all ijj E Sy'- they are exactly the subspaces 
V with finite average j^^^ T{ip) dili due to the relation HanV^H^ < ||an||F = J2k ll^n^fcP 
for any state ip and orthonormal basis (pk of V. This means that the presence of a 
single state ip E Sy with T{ip) = oo implies the divergence of the average j^^^ T{ip) dip. 

The reason for this is that the set of states if) E Sy with finite t{iI)) has the 
structure of the unit sphere of a subspace of V . This follows from the inequality 
\\an{ai) + P(t))f < 2|an|a„^f + 2|/3|2||a„0f, which gives r(a^ + /30) < 2|a|V(^) + 
2|/3pr(0). Therefore, unless this set exhausts all of Sy, it has null Lebesgue measure 
with respect to the uniform measure dip. Thus, T{ip) diverges dip-a-.e. on Sy, which 
causes the divergence of the average j^^^ T{ip) dip. 

The following theorem characterizes the subspaces V with a finite averaged ex- 
pected ^-return time and gives a formula for this average. It can be considered as 
the extension to subspaces of the results given in [151 Theorem 2]. A key ingredient 
will be the determinant det T of an operator T on V, that is, the determinant of any 
matrix representation of T. 

Theorem 3.2. Consider a unitary step U and a finite- dimensional subspace V with 
spectral measure ^{dX), Schur function f{z) and first V -return generating function 
a(z) = zf\z). Then, the following statements are equivalent: 

(1) All the states ofV are V -recurrent with a finite expected V -return time. 

(2) All the states of V are recurrent with a finite expected return time. 

(3) fi{dX) is a sum of finitely many mass points. 

(4) f{z) is rational inner. 

(5) det f{z) is rational inner. 

Under any of these conditions, the average of the expected V -return time is 

with K a positive integer that can be computed equivalently as 

K = J2^MEkV) = ^rank^({Afe}) = deg det a(e^^), (9) 

k k 

where Xk are the mass points of n{dX), = E{{Xk}) are the orthogonal projectors 
onto the corresponding eigenspaces of U = J XE{dX) and deg det a{e^^) is the degree 
o/deta(e*^): — )■ S^, i.e. its winding number, which coincides with the number 
of the zeros of det a{z) inside the unit disk, counting multiplicity. 

Proof. We know that (3) -v^ (4). On the other hand, the polarization identity implies 
that n{dX) is finitely supported if and only if the spectral measure = {ip\E{dX)ip) 
of any state ip E V is finitely supported too, so (2) (3) by [HI Theorem 2]. 
Hence, to prove the equivalences we only need to argue that (1) =^ (5) =^ (4) =^ (1). 
The formula T^ip) dip = ^.^ y will be established as a key step of the proof that 

(4)^(1). 
(1) ^ (5) 
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The recurrence of V means that f{z) is inner, i.e. it is unitary a.e. on the unit 
circle. Therefore, det f{z) is analytic on the unit disk with unimodular boundary 
values, which implies that it is a scalar inner Schur function. 

The finiteness of the average for the expected V^-return time can be expressed by 
requiring that /(e*^) belongs to the space iJ^/^ of operator-valued functions h{9) = 
EZ-oohnC''^' such that \\h\\l,,, = E^=-ool^lll^n||| < oo. Siuce /(e*^) G L^, it 
follows that the matrix element {4>\f{z)il)) lies in the scalar H^^"^ fl L°° space for any 
states (jy^ifj E V . Bearing in mind that H^/'^ fl is an algebra, we conclude that 
det/(e^^) G H^/^. 

Hence, det/(2;) is a scalar inner Schur function in H^^^. On the other hand, the 
one-dimensional analogue of this theorem, already proved in [151 Theorem 2], implies 
that a scalar inner Schur function is in H^^^ if and only if it is rational. Therefore, 
det f{z) is rational inner. 

(5) ^ (4) I 

As any rational inner Schur function is a finite Blaschke product, 

d 

det/(z) = Cn^(^)' h^{^) = ^—^. kl<l, ICI = 1, 

with finitely many zeros Wi, . . . , inside the unit disk. 

We may factorize these zeros in the operator Schur function f{z) by using the 
rational inner Blaschke-Potapov factors [H] 

B^,^{z) = Iv + (b^z) - E Sv 

From Potapov's factorization, for some unit vectors (pj G ker f{wj), 

f{z) = -B^i,0i(2) ■ ■ ■ 9{z), 

with g{z) analytic on the unit disk. Since g{z) must be unitary on the unit circle 
too, it is an inner Schur function. 

Bearing in mind that det Bu,^^{z) = bu,{z), we find that det g{z) = C is a unimodular 
constant for \z\ < 1. Consequently, g{z) has an analytic inverse on the unit disk 
and, since the values are unitary on the boundary, not only < 1, but also 

||fl'(^)~"'^|| < 1 for l^;] < 1. This implies that = 1 for any E Sy and 

|z| < 1. Hence the polarization identity and the maximum principle show that g{z)ip 
is constant on the unit disk. We conclude that g{z) is a constant unitary, therefore 
f{z) is rational inner. 

(4) (1) I 

It suffices to derive the formula for j^^^ t(^/') dtp since (1) is equivalent to the finite- 
ness of this average. Any rational inner Schur function f{z) is analytic in a neighbor- 
hood of the closed unit disk. Therefore, the restriction of a,(z) = zf\z) to the unit 
circle is a unitary function a(e*^) that is real analytic in 9. Hence, the expression of 
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the average of r (■?/') in terms of radial limits to the unit circle becomes 

/ r(^) = K = r Tr (a(e^'')t9,a(e^^)) ^. 

Jsv Jo 27r2 

The extension of Rellich's theorem [221 [13 to real analytic families of normal op- 
erators ensures the existence of an orthonormal basis of V consisting of real analytic 
eigenvectors 0^(6*) of a(e*^) with real analytic eigenvalues C,k{(^)- 

Differentiating a{e'') = Efc6(^)l0fc(^))(0fc(^)|, we get 

a'^dea = '^^kde^k\h){4>k\ + ^ {^j^k{4>j\d04>k) + {d0^j\4>k)) \4>j){4>k\, j^xq) 
k j^k 

where we used Re(0fc(6')|900fe(6')) = (since = 1). This implies 

Tr {a{e'ydea{e'')) = J2^9gU0) = det a(e*^) de det a{e''). 

k 

Therefore, 

K = / = deg det a(e'^) 

is the winding number of deta(e*^): — t- S^, which is finite because it coincides 
with the number of zeros of the finite Blaschke product det 0(2;). The alternative 
expressions given in for the positive integer K are proved in Theorem IA.2I of 
Appendix El D 

The previous theorem states that the presence of a finite-dimensional subspace V 
whose states return to in a finite expected time is equivalent to the existence of 
eigenvectors for the unitary step U. However, as mentioned previously, when such 
states do not cover all of Sy, no singular subspace is required for U. For instance, 
in |0) returns to V = span{|0), |1)} in one step with the unitary shift U\x) = 

|a; + 1), a; G Z. 

The relation ( llUp sheds light on the origin of the non-integer values of the geomet- 
rical phases T{ip) and the reason for the cancellation of the non-integer contributions 
when averaging over the states tp & Sy- The identity 

(^|r/)(0|V;)#= V0,r/GH, 
Sv dim V 

shows that such an average cancels the off-diagonal terms in flTU]) and transforms each 
projector |0fc(^))(0fc(6')| into the constant 1/ dimy. This eliminates the dependence 
on 0fc(^) so that the only contributions to 't{'4') dip come from the winding numbers 
of the eigenvalues ^k{d)i which are integers. However, if we do not perform the 
average, the transport of the basis 0fc(^) along the closed path gives off-diagonal 
contributions, and also modifies the diagonal ones due to the ^-dependence of the 
projectors |0fe(^))(0fc(6')|, causing the emergence of non-integer values. Nevertheless, 
T{'ip) is independent of the choice of basis functions 0fc(6') because the definition of 
the expected V^-return time does not require such a choice. In other words, while the 
expected l^-return time t{iI)) is a non-integer geometric property of the loop in Sy 
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defined by ip{9) = a{e^^)ip, the average over Sy of these geometrical phases leads to 
a true topological invariant, the winding number of deta(e*^), which is therefore an 
integer. 

Concerning the meaning of the requirement that fi{dX) is finitely supported, re- 
member that /x((iA) = PE{dX)P, so the mass points of fi{dX) must also be mass 
points of E{dX), which means that they are eigenvalues of U. If fi{dX) is a sum of 
finitely many mass points A^, then P = P J E[dX)P = P^^E^P, which is equiva- 
lent to stating that V C ®kEkV. Indeed, if ijj e V and P = PQP with Q = Y.k ^k, 
then ip = PQip. So Qip = ip + (p with G V-^ . \\Q\\ = 1 implies = 0. Conversely, 
if every ip & V may be written ip = ^j^Ekipk with ipk € V, then Qip = ip for all 
ip E V. Hence, PQP = P. In summary, since Ek is the orthogonal projector onto 
the eigenspace associated with A^, we find that the spectral measure of V is finitely 
supported if and only if is a subspace of a finite sum of eigenspaces of U. 

The [/-invariant subspace "Ho = ®kEkV whose dimension gives K has a remarkable 
meaning. Consider an arbitrary [/-invariant subspace "H. The restriction If = U \ 
?^ of [/ to ?^ is a unitary with spectral measure E{dX) = E{dX) \ H and thus 
with eigenspaces Vk = Ekl-i. Hence, ii V <Z l-i, then EkV C V^, which implies 
that I-Lq C T-L. This identifies T-Lq as the minimal [/-invariant subspace in which 

V lives, and allows us to state that the average of the expected V -return time is the 
ratio (\im.{®kEkV) / dim V between the dimension of the minimal U -invariant subspace 
containing V and the dimension of V . 

In the case dimH < oo the integer K can also be calculated as dimH — where 

V is the number of linearly independent eigenvectors of U lying in V-^ (see Appendix 

m- 

An additional consequence of the previous theorem is an identity that relates the 
degree of the determinant of a D-dimensional rational inner Schur function f[z) and 
the masses of the corresponding measure fj,{dX). Since a{z) = zf\z) we have that 
det a{z) = z^det f{z), thus 

degdet/(e*^) = rank /x({A,}) - /}, 

k 

where A^ are the mass points of iJ,{dX). 

4. Applications: site recurrence 

The most natural example of subspace recurrence is what we can call site recurrence 
of a quantum walk. Consider a quantum walk on a spatial lattice A with a set of 
internal degrees of freedom E, i.e. the basis states (l)x,a are labelled by the site index 
X G A running over the lattice and an additional index a G S that refers to the 
possible basis states for each site. If we consider the evolution starting from a state 
ip = Ca(px,a localized at a given site x, we can ask for the return probability 

to the original site. In other words, we deal with the notion of Kj-recurrence for the 
subspace = span{02;,a}«gs- 

In this section we will present explicit analytical results for the site recurrence of 
any coined quantum walk on the line, the half-line and a finite one-dimensional lattice 
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too. Some numerical calculations (given later) will show the effect of the dimension 
and the geometry of the lattice on site recurrence. 



4.1. Site recurrence in one dimension. Coined walks on the line live on the 
Hilbert space ® C^, whose basis states we write as |x, a), x G Z, a G {t, i}- The 
step operator U = SC is factorized into a shift S\x,^) = |x + 1, t), S\x,l) = \x — l,i) 
and a coin operation C acting at each site subspace = span{|x, t), \x, I)} separately 
as a two-dimensional unitary operator. Bearing in mind the phase factor freedom of 
the basis states, we can assume without loss that the matrix representation in the 
basis t), l^; i)} of the coin operator at site x has the form 

a = ((x,a|C|a;,/3))„«6|t ,| = ( , " 

The convenience of using even indices for the coefficients of Cx will become clear later 
on. 

Actually, we will assume the strict inequality 172x1 < 1- Otherwise the quantum 
walk would decouple trivially into independent ones hving on orthogonal subspaces 
generated by basis vectors that could be analyzed independently. 

As we will show, the study of site recurrence for one-dimensional quantum walks 
is greatly aided by the theory of orthogonal polynomials on the unit circle and CMV 
matrices; (for a very quick and modern introduction to both topics see [28] )• Indeed, 
they provide efficient tools for the analytic computation of the 2x2 Schur functions 
that encode the site recurrence. Note that the polynomial approach becomes simpler 
for quantum walks on the half-line or a finite one- dimensional lattice than for the case 
of the whole line [H [9l [10]. In the former case, the Schur function is scalar- valued, 
while in the latter case, it is matrix-valued. 

Our first step will be the study of site recurrence for coined walks on the half-line 
with reflecting boundary conditions at the origin. They are defined as in the case of 
the whole line simply changing the lattice Z by = {0,1,2,...} and the shift at site 
by 5*10, 1) = |0, t)- The case of a finite lattice = {0, 1, . . . , A^ — 1} with reflecting 
boundary conditions at the edges simply requires an additional modification of the 
shift at the right edge by S\N - 1, t) = |A^ - 1, i). 



4.1.1. Coined walks on the half-line. According to the previous paragraph, the tran- 
sition matrix of a coined walk on the half-line has the form 



C = i{ej\Uek))j, 



( 7o 


Po 
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Po 
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-74 
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where we assume the ordering e2x = 1^;, t)) ^2x+i = for x G Z+. This matrix is 
the special case 72^+1 = 0, /c G Z+, of the more general unitary 

( lo Po 

Poll -loll P1I2 P1P2 

PoPi -loPi -111 2 -71/52 

P273 -7273 P374 P'iPi 

P2/53 -72/53 -7374 -73/54 

P475 -7475 

P4/95 -74/55 

known as a CMV matrix with Verblunsky parameters (7^) [TTl [T2ll26| [28 | [29 | IM] . 

CMV matrices are closely related to the orthogonal Laurent polynomials on the 
unit circle, X}~{z), which are defined by performing the Gram-Schmidt procedure on 
the set {1, 2, 2"^ 2^, 2;"^, , . . .} m. for a given measure fi{d\) on S^. They satisfy 
the recurrence 

X{z) C = zX{z), X{z) = {Xo{z) Xi{z) X2{z) . . .) , Xo(^) = 1. 

The CMV shape of the transition matrix provides a direct connection between 
coined walks on the half-line and orthogonal Laurent polynomials on the unit circle 
with null odd Verblunsky parameters. For more details, see [H]. 

The Verblunsky parameters of a CMV matrix have another important mathemati- 
cal role: they characterize the Schur function f{z) of the orthogonality measure fi{dX) 
via the Schur algorithm 

/oW = /(^), 

/.+i(^) = - /'^j;7v 7. = /.(o), k>o, 

which generates a sequence /^(z) of Schur functions called the iterates of j{z) [23 | [2 ^ 
[IHI, |26l [28] . Geronimus' theorem [26] implies that the Verblunsky parameters (7^) in 
the CMV matrix are equal to the Schur parameters (7^) in Schur's algorithm. Note 
that the /c-th iterate jh{z) has Schur parameters (7fc, 7fc+i, 7fc+2, • • • )• 

Coined walks on the half-line correspond to Schur parameters (70, 0, 72, 0, 74, 0, ... ). 
From the Schur algorithm, we see that the absence of odd terms is equivalent to the 
Schur function /(z) being an even function of z. 

The orthogonality measure of a CMV matrix has an infinite support iff there are 
infinitely many Schur parameters. If the measure is supported on a finite number of 
points, say A^, the corresponding Schur function will have the same finite number of 
Schur parameters (70, 71, ... , 7Ar-i), all of them in the open unit disk except the last 
one, which lies on the unit circle. The reason for this is that the Schur function is a 
Blaschke product of degree TV — 1, i.e. a rational inner function, and thus the Schur 
algorithm stops after — 1 steps because it reaches the unimodular constant function 
7Ar_i. Therefore, an inner Schur function is rational exactly when it has finitely many 
Schur parameters. 



\lk\ < 1, 

Pk = V^-hk 
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Together with the standard iterates, the inverse iterates are key in the study of 
recurrence. If the Schur function f{z) has Schur parameters (70, 71, . . . ), the rational 
inner Schur function f''{z) with Schur parameters (— 7^, — 7fe_i, . . . , — 7o, 1) is called 
the k-th inverse iterate of f{z). This terminology was introduced in [19], where the 
author uses the notation bk-i{z) for the function we denote f'^{z) for convenience. 

An important consequence of the above results is that the step operator [/ of a 
coined walk on the half-line is unitarily equivalent to the multiplication operator 
on given by h{z) 1— )■ zh{z). The equivalence is estabhshed by H- Xk{z), so 



(ej|f/"efc) = / Xj(A) A'*Xfc(A) /i(rfA) for n e Z. On the other hand, the spectral 
decomposition U = J \E{dX) yields (ej|f/"efc) = / X^{ej\E{dX)ek). This identifies 
the complex- valued measure 



in terms of the orthogonal Laurent polynomials and the corresponding orthogo- 
nality measure. Hence, the spectral measure fi^{dX) = {ip\E{dX)ip) of any state 
^ = I]fc>oCfc^fc l^-^idX) = \h{X)\^fi{dX) with h{z) = Y.k>o'^kXk{z) e L^, 

so it is absolutely continuous with respect to ^{dX). In particular, the orthogonality 
measure fi{dX) coincides with the spectral measure of the first basis state Cq = |0,t)- 
Concerning site recurrence, identifying the spectral measure ix^{dX) of the site 
subspace = span{e2x, ^2x+i\ with its matrix representation in the basis {622;, ^2x+i}^ 
we can write 



where fk{z) and f^{z) are the iterates and inverse iterates of the Schur function 
f\z) corresponding to the orthogonahty measure ^{dX). More precisely, the Schur 
parameters of f2x+i{z) and p^~^{z) are 

(0,72x+2,0,72x+4, •••) and (0, -72^-2, 0, -72x-4, -7o, 1), 

respectively. In particular, f^^~^{z) is a finite Blaschke product of degree 2x. 

To obtain the site return probability operator we must simply take into account 
that f^{z) is an inner function. Then, the previous matrix Schur function gives the 
Vir-return probability matrix 



{ej\E{dX)ek) = X,{X)Xk{X)fi{dX) 



H^{dX) = {{ej\E{dX)ek))j,ke{2x,2x+i} 




(12) 




(13) 




f |72xP + pL|I/2x+iP P2x72x(l - ||/2x+lin 
VP2x72x(1 - \\f2x+lf) pL + l72xP||/2x+lf 

/2 — (1 — ||/2x+iinr2;, = [ I 




) 
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where I2 is the 2x2 identity matrix and || ■ || is the norm with respect to Lebesgue 
measure on the unit circle. 

The matrix has eigenvectors (72., p2xY ^i^nd (p2x, —l2x)^ with eigenvalues and 
1 respectively. Therefore, the maximum and minimum return probabilities to the site 
X are 1 and ||/2x'+i|p5 attained respectively at the qubits 

V'i^^ =72x|x,t) + P2x|a;,;), = p2x\xA) -'^2x\x.\). (14) 

The average of the return probability to the site x is therefore 

I I f ■ _ I —I— II 10.. I -1 



i 



TiR., I+II/2X+1I 



Sv, dim K 2 

From the above results, we find that any site of a coined walk on the half-line has 
at least one site-recurrent qubit, i.e. a qubit that returns to the site with probability 
one. The orthogonality measure fi{dX) of a CMV transition matrix like f|TT]) can be 
absolutely continuous, from which it follows that the spectral measure fi^{dX) of any 
state is also absolutely continuous. Therefore, coined walks on the half-line associated 
with absolutely continuous orthogonality measures p{dX) are examples of quantum 
walks without a singular subspace for the unitary step U, but having a l^-recurrent 
qubit for any site subspace V^. 

Furthermore, if a site is recurrent then all the states of the Hilbert space must 
be recurrent, which implies that all the sites are simultaneously recurrent. Indeed, 
the previous results show that is a recurrent subspace if and only if f2x+i{z) is 
inner. Since a Schur function f{z) and its iterates fk{z) are simultaneously inner or 
non-inner, we conclude that the presence of a recurrent site implies that the orthog- 
onality measure fi{dX) is purely singular. Then, not only does every qubit become 
site recurrent, but every state becomes recurrent as well since its spectral measure is 
absolutely continuous with respect to fi{dX). 

A remarkable fact is that the return probability matrix for the site x only 
depends on the Schur parameters (72X, 72x+2, 72a;+4, • • • )• This means that the site 
return probability of a qubit at a site x does not depend on the coins Cy, y < x, of 
the previous sites. This is the site extension of a result already found in [13 Section 
6.1] for the state recurrence of coined walks on the half-line. 

Concerning the expected return time to the site x, we should distinguish two cases 
depending on whether the measure fi{dX) is purely singular or not. If not, the expected 
V"a;-return time makes sense (i.e. it is not necessarily infinite) only for ip^x^ because it is 
the unique site-recurrent qubit at site x. We find that r(?/'i^'') = 1 + 5r(^/'i"'^^) = 2a; + 1 
is an integer because 

5r(V;«)=lim ^(72. P2x)deUre'')Ure'y(^''' 



dep-\e^') dB_ ^ 

i.e. 5t{'iI)^'^) is the winding number of the finite Blaschke product /^^~^(e*^) of degree 
2x. 
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If fi(dX) is purely singular, then all the sites are recurrent and the expected return 
time to the site x makes sense for any qubit ipx = C(\x,t) + (3\x,i). Just as in the 
previous case we obtain 



T{ipx) = 1 + |a72x + /3p2xpdeg/ 



2x-l 



+ lim \ap2a 

r— 



/372xpTr(/2x+l), 



where, for any Schur function g{z) = J2n>o^-n-^^ ^ (^j l); use the notation 



Trig) 



2n 



de 



E 

n>l 



(15) 



When g{z) is rational inner, obviously \im.j.^i- Tr{g) = degg. On the other hand, 
if g{z) is inner but not rational - i.e. it is inner with an infinite sequence of Schur 
parameters - it is proved in [151 Theorem 2] that lim^^i- Tr{g) = oo. This last one is 
the situation for the Schur function f{z) of fi{dX) and its iterates. Hence, we finally 
get 

2x + 1, if ipx e spa.n{ilJx''}, 
oo, otherwise. 



As a consequence, j^^^ t(^) dtp = oo, which is in agreement with the fact that the 
spectral measure fi^{dX) in f[T^ is never finitely supported because the orthogonality 
measure p{dX) has an infinite support. 

These results prove that for any site x, coined walks on the half-line always have a 
\4-recurrent state with a finite expected ^.-return time, independently of the presence 
or absence of a singular subspace for the unitary step U. 

4.1.2. Coined walks on a finite lattice. The ordering cq = |0,t),ei = |0, |), . . . , e2Ar_2 = 
I — 1, 'I"), e2N-i = \N — 1, D of the basis states yields the following transition matrix 
for a coined walk on the finite lattice Z^r, 



/7o 


Po 



c 



Po 

72 
-70 


P2 



P2 




-72 



V 



72Af-4 




P2N-A 



P2N-A 




-12N 



l2N~2 
2 
P2N-2 



P2N~2 



-72Ar-2 



This transition matrix is a finite CMV matrix with Verblunsky parameters 

(70,0,72,0,...,72Ar_2,l). 

Unlike the case of a semi-infinite CMV matrix, this one generates only a finite segment 
of 2N orthogonal Laurent polynomials. The orthogonality measure fi{dX) is again the 
spectral measure of the first basis state cq, but now it is supported on 2A^ points. 
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Consequently, for a coined walk on a finite one- dimensional lattice, every subspace is 
recurrent with a finite expected return time for any state. 

Concerning site recurrence, it remains to calculate the expected return time for 
any qubit. For this purpose we need the Schur function f^^iz) of a site x, which is 
the same one already given for the half-line, but taking into account that the iterate 
f2x+i{.z) of the Schur function f{z) for fi{dX) now has a finite segment of Schur 
parameters (0, 72X+25 0, 72^+45 • • • , 0, 72Ar-2, 1), and thus it is a finite Blaschke product 
of degree 2(A^ — 1 — x). We conclude that the expected V^-return time for a qubit 
ipx = Ci\x, t) + /3\x, I) is given by the quadratic form 

a 



t{^Px) = {a P)'rxip 



with Tx = h + Stx and 

Srx= r deUe'')fx{e''''- 







\l2x? deg + pI^ deg f2x+i P2x72x(deg p""'^ - deg f2x+i] 
P2x72x(deg /^^-^ - deg f2x+i) pL deg P""'^ + |72xP deg f2x+i^ 

= degp^-'h + (deg/2.+i - deg p^-')T,. 

The extreme values of the expected return time for the site x, attained at the 
qubits vi^^ and ijj^^ given in f lT^ . are the integers 1 + deg p^~^ = 2a: + 1 and 
1 -|- deg /2X+1 = 2(iV — 1 — a;) -|- 1, which do not depend on the coins but only on the 
number of sites to the left and the right of the site x respectively. However, there 
exist qubits at site x whose expected return time to such a site is any value between 
the above extreme integers, which shows the non-integer nature of the Aharonov- 
Anandan phase giving the expected V^-return time. Despite this fact, the average of 
t{iPx) is, not only an integer. 



i 



TtTx (l + deg/^--i) + (l + deg/, 



r(^) = -^—^ = , , , ^^^j2x+i) ^ ^ 



Sv, dim Vx 

but also independent of the site x and the coins of the walk: it is always equal to the 
number of sites of the lattice. 

The above average can be calculated also from the degree of 

detd,(^) = det^/](^) = -z^P^{z)7^i{z), 

which gives 

f ( i I deg det 2 + deg /^^"^ + deg /z^+i , , 

4.1.3. Coined walks on the line. The transition matrix with respect to the basis or- 
dered as e2x = \x, t), G2x+i = \x,D, X & Z, is now a doubly-infinite CMV matrix with 
Verblunsky parameters (. . . , 0, 7_2, 0, 70, 0, 72, 0, ... ) [SI 13 flU] . 

As in the previous cases, the analysis of site recurrence for the whole line comes 
from the Schur function of a site. Appendix |B] shows that this Schur function can be 
read from the case of the half-line by extending to a doubly-infinite sequence of Schur 
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parameters (7fc)fcgz the definitions of fk{z) and /'^(z) as tlie Scfiur functions with 
infinitely many Schur parameters (7^, 7^+1, 7^+2, • • • ) and (-7^, -7^-1, -7^-2^ • • • )' 
respectively. In other words, the expression for /^(z) given in f|T3|) for the half-line 
also holds for the whole line with f2x+i{z) and /^^"^(z) corresponding to the Schur 
parameters (0, 'J2X+2, 0, 72^+4, 0, . . . ) and (0, -%^_2, 0, • • • ) respectively. 

The main difference compared to the half-line case is that p^~^{z) is not a rational 
inner function on the line. Therefore, the return probability matrix to the site x takes 
the form 

R _ (\l2x\V~'r + pUf2x+ir P2.72x(||/^^^i^-||/2x.+lin 

' ~ \P2x%xi\\p''~'r - ii/2.+iin pUf^-'r + 172.P11/2.+1IP 

= ll/'^"'lP/2 + (||/2.+il|-||/'^-^f)r.. 

The extreme values of the V^-return probability are ||/2x+i|P and they are 

again reached at the qubits ipx''' given in (fT^ . The average over of the return 
probability to the site x is 



Site-recurrent qubits appear if and only if p^~^{z) or f2x+i{z) are inner for some 
site X, and a site x is recurrent exactly when p^~^{z) and f2x+i{z) are both inner. 
However, the fact that fk+j{z) and f^~^{z) are the j-th iterates of fk{z) and f^{z) 
respectively implies that the Schur functions /2x+i(^) are simultaneously inner or 
non-inner for all x, and analogously for p^~^{z). Hence, concerning site recurrence, 
all the sites exhibit the same characteristics: either none of them has a site-recurrent 
qubit; or all of them have a single site-recurrent qubit; or all of them are recurrent 
sites. 

Further, as on the half-line, the presence of a recurrent site implies not only the re- 
currence of any other site, but the recurrence of any state, and hence of any subspace. 
The reason for this is that the spectral matrix measure fi{dX) of the subspace spanned 
by the basis states |0, t), | — 1, i) generates the measure of any other state through an 
expression of the form h{Xy fi{dX)h{X) for some vector function h{X) G L^, while the 
matrix Schur function of fj.{dX) is ITU] 



fo{z] 

r\z) 



Obviously f2x+i{z) and fo{z) are simultaneously inner or non-inner, as are p^~^{z) 
and f~^{z). Hence, the presence of recurrent sites is equivalent to fo{z) and f~^{z) 
both being inner, which at the same time is equivalent to f{z) being inner. This 
means that the matrix measure ii{dX) - and thus the spectral measure of any state 
- is purely singular, implying the recurrence of any state. 

A special situation holds when ||/^^~^|| = ||/2x+i||; a condition that characterizes 
the sites x with an effectively scalar site return probability matrix = ||/2x+i|P-^2; 
i.e. with the same site return probability for any qubit. This holds for instance if the 
Schur parameters are symmetric (72x-2fc = e*'^''"'"^'^'^^722;+2fc) or conjugated {'y2x-2k = 
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^i{r]+2kui)^^^^^^^ with respect to the site x, up to a hnear phase. In the former case, 
/2^-i(^) = -e-''^J^{e-''^z), and in the latter, P'^-^z) = -e-'"^ f2cc+i{e-"^z), as 
follows from the Schur algorithm. In particular, for a constant coin on the line, all 
the qubits of all the sites have the same site return probability since f2x+i{z) and 
P^~^{z) = —f2x+i{z) are independent of x. 

Given a site-recurrent qubit ipx = f) + f3\x,l), its expected return time to the 
site X is given by 

I2 + dTx,r, 



p2xl2x{Tr{P'' ^) - Tr{f2x+l)) \ 

plxrr{P^^') + \l2x?rrU2x+,)) 

PIx_ -P2xl2x\ 
P2xl2x l72xP y ■ 

where was defined in f lTSj) above. The matrices and are positive semidefinite 
and their kernels - spanned by ■j/^i^'' and -j/^i^^ respectively - are orthogonal to each 
other. Note that rr(-) > 0, and lim^.^!- Tr{fk) = 00 whenever fk{z) is inner because 
it has infinitely many Schur parameters, with a similar result for f^{z). As a conse- 
quence, T^ipx) = 00 for any site-recurrent qubit ipx- Thus, in contrast to the case of 
the half-line, the expected return time to a site is infinite for any site-recurrent qubit 
on the whole line. 

We have seen that the analysis of site recurrence for any coined walk on a one- 
dimensional lattice can be reduced to the study of scalar Schur functions, despite the 
fact that subspace recurrence is characterized by matrix Schur functions. This is due 
to the form of the site Schur function, which is given in terms of 'right' and 'left' 
scalar Schur functions characterized by the Schur parameters of the coins located to 
the right and left of the site under consideration, respectively. Moreover, a reflecting 
boundary condition at a left or right edge implies that the corresponding scalar Schur 
function has a finite number of Schur parameters, and thus is rational inner. This 
results in the independence of the site return probability with respect to the coins on 
the side of the site where the reflecting boundary condition appears. In any case, the 
return probability matrix to a site only depends on the coin at that site and the 
norm of the right and left Schur functions. 

These norms can be computed in some cases, providing a completely explicit ex- 
pression for the site return probability matrix. For instance, this is the case for a 
constant coin 

C=(^^ 7) , P = v^T^, (16) 



Stx 



hm (a /3) Tx,r , Tx,r 

Ztt 

72xP TriP""-^) + rrif2x+l. 
P2xl2x{'rr{P''~^) - Tr{f2x+l)) 

^^,(/2-l)r- + r,,(/2.+l)r., 

\l2x? P2xl2x 
P2xl2x pIx 



28 



J. BOURGAIN, F.A. GRUNBAUM, L. VELAZQUEZ. J. WILKENING 



which yields [IHl Hn] 

||/2.+if = 11/^^-1' = ^ {p|7l + (1 - 2/) arcsin I7I} 

7r|7p 

for any site x, assuming no reflecting boundary conditions at the right or left edge, 
respectively; (otherwise the related norm is one). 

Finally, from the Taylor expansion around the origin of /2x+i(^) and we 
can obtain the Taylor coefficients of the site Schur function f^iz) = ^„>i al^^z"'~^, 
which provide the n-step first return amplitude matrix an^x to the site x. For in- 
stance, the constant coin (fT6!) on the line yields f2x+i{z) = J2n>i'^nz"' and p^~^^( 

-E„>iCn^" with danu] 



z 



ci=7; C2„ = 0, C2n+i = — — — , c=l-27 , n>l, 

27(72 + 1) 

and Pn{z) the Legendre polynomials. Hence, 

0'x,2n-l = 0, Ctx,2n = C^n"^, ^ > 1, 

Therefore, given two qubits ipx-, 4'x at site x, the n-step amplitude for the transition 
4'x without hitting the site x at intermediate steps is dn{4>x\'^'4'x) ■ The time 
and qubit dependence of these amplitudes enter through the scalar and matrix factors 
dn and T, respectively. T encodes the ratio between the n-step amplitudes (without 
hitting the site earlier) for different pairs of qubits, which are independent of n and 
the site. The time dependence, isolated in the scalar factor dn, shows an n~^/^ decay 
as n oo [IS [TU] . 

The comparison between state and site recurrence for coined walks leads to inter- 
esting phenomena seen previously in Example 12.31 We can compute the state return 
probability of a qubit following steps similar to those in Example 12.31 Indeed, us- 
ing the matrix Caratheodory function Fx{z) corresponding to the site x, obtained 
from the site Schur function fx{z), we can arrive at the scalar Caratheodory function 
(a0) Fx{z) (^) of any qubit ipx = t) + f^\x,i)- The Lebesgue norm of the 
related scalar Schur function yields the return probability to the same qubit. 

The upper plot in Figure [2] shows the state and site return probabilities for a = cost 
and /3 = sint in the case of a constant coin on the half-line given by 72^ = 
We see that this kind of coined walk exhibits the same paradoxical behavior already 
found in the toy model of Example 12.31 since the state return probability is higher 
than the site return probability for some values of t. By contrast, this effect is absent 
in the lower plot of Figure [21 which represents the analogous situation on the line for 
a = cos t and /3 = i sin t. 

Indeed, for any constant coin 72^ = 7 on the line, the above procedure for the 
computation of the state return probability for a qubit ipx = C(\x,t) + f3\x,i) yields 
Ra,i3 = \\f{z,0)f{z,c)\\'^, where fi^z.-^o) is the Schur function with Schur parameters 
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0.44 



0.4 J L 

0.25 0.5 0.75 1 

Figure 2. State versus site recurrence on the half-line and 

the line. Given the qubits ipx(t) = cost|x,t) + sint|x,|), the upper 
figure represents as a function of t the state (blue) and site (red) return 
probabilities for a constant coin 72a; = -y/S/S on the half-line. The lower 
figure is a similar plot for the same coin on the line, but for the qubits 
ipxit) = cost |x, t) + isint \x,l). 



(70,0,7,0,7,0,...) and c = 7 — (2zp/7) Im(a/37). On the other hand, our previ- 
ous results imply that the corresponding site return probability has the same value 
{ipxlRx'ipx) = ||/(^)0)|p for any qubit. Thus, for any constant coin on the line, the 
state return probability of a qubit is never higher than its site return probability. 
Nevertheless, this can change if we allow distinct coins on the line or if we consider 
absorbing subspaces not attached to a site. 

4.2. Site recurrence in two dimensions. In this section we show some numeri- 
cal results and display a few figures illustrating the behavior of return probabilities 
starting from different states in a few types of two-dimensional coined walks. 

We first consider the usual two-dimensional doubly-infinite square lattice k = 1? 
with four internal degrees of freedom for each site (east, north, west and south). 
The Hilbert state space £^(Z^) ® is spanned by a basis of the form |x, — )■), |x, t), 
|x, ^), |x, J,), X G . The unitary step U = SC factorizes again into a shift S\x, — )■) 
= \x + ei, -)■), t) = \x + 62,1), 'S'lx, ^) = \x - ei, ^), S\x,l) = \x - 62,1), with 
61 = (1,0) and 62 = (0, 1), and a coin operation C that acts independently on each 
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site subspace 14 = span{|a:, — )■), |x, f), |a:, |)} as a four-dimensional unitary- 

operator. 

The walks we consider are the well known Grover walk and the Fourier walk, defined 



respectively by the following 


constant coins in the site basis, 






/-I 1 


1 1\ 




/l 1 1 


1 \ 


^ 2 


1 -1 
1 1 


1 1 
-1 1 




1 I -1 
1 -1 1 


—i 
-1 




U 1 


1 -V 




\l -i -1 


' J 



The method we use to produce these results is a classical one. It goes by the 
name of the "transfer matrix method," and our best references for this material are 
the articles |6l [7] and the classical book [30]. The main point is that once the four- 
dimensional unitary coin is specified, one has a way of obtaining the four-dimensional 
return amplitudes for any site. The methods of this paper then allow one to 
compute many other quantities of interest. 

Figure [3] compares the state and subspace return probabilities computed on a curve 
of qubits lying in the same site. Each figure represents the return probabilities to 
several nested absorbing subspaces of the site subspace to reveal the possible exotic 
behavior of the return probability. Although the Grover example shows a classical 
monotone behavior, it conceals an interesting effect because the return probability 
to the four-dimensional site subspace coincides with the return probability to one 
of its three-dimensional subspaces. Furthermore, the two examples of the Fourier 
walk show that the return probability is not a monotone function with respect to the 
inclusion relation among absorbing subspaces. 

To see the effect of the geometry of the lattice in the return probabilities we also 
present in Figure H] a similar comparison for a couple of coined walks in a two- 
dimensional hexagonal lattice already used by other authors for different purposes 
|20] . In this case each site has three internal degrees of freedom that control a con- 
ditional shift S to the three nearest neighbors analogously to the case of the square 
lattice. A three-dimensional unitary operator at each site defines a coin operation C 
giving the unitary step U = SC. The examples represented in Figure H] correspond to 
the constant coins given in the site basis by the three-dimensional versions of Grover 
and Fourier, 

, /-I 2 2 \ 

These figures show that the exotic dependence of the return probability on the absorb- 
ing subspace holds for two-dimensional coined walks with a constant coin for many 
lattice geometries. Indeed, the plots of Figure Hlshow more dramatic phenomena than 
the previous ones. 

Each of these two-dimensional coined walks has a constant coin and no boundary 
conditions; thus, the site return probability operator R = is independent of the 
site X. The computation of its eigenvalues reveals a curious doubling effect: there is 
always an eigenvalue with multiplicity two. 
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Figure 3. State versus site recurrence on a square lattice. 

For a choice of qubits (pl, (pl, 4>l at site x, the figures sliow as a 
function of t the state (blue) and site (black) return probability of 
i'xit) = cost 4)1 + sin t (pl, together with the return probability to the 2- 
dimensional subspace span{0^, 0^} (green) and the 3-dimensional sub- 
space span{0^, 0^, 0^} (red). The upper figure represents the case of the 
Grover coin with (pl = \x, -)■), (pl = \x, (pl = -^{\x , "\-) + i\x , i)) . 

Curiously, the change (pl = -^{\x,^) + \x, D) yields a 3-dimensional 
subspace with the same return probability plot as the whole site. The 
other two figures correspond to the Fourier coin Cp with the choice 
(1^1 = -^), 0x = k.t), 0x = 73(k'^) + k'i)) (center) and 
<Pl = + (bottom). 
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— state / 




— dim/ ^ — \ 

— site ^^^v^ 
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0.25 0.5 0.75 1 



Figure 4. State versus site recurrence on a hexagonal lattice. 

Different return probabilities as a function of t for the qubits given in the 
site basis by V^^l^) = 75 ^ ^' + « sint (0, 1, 0). Together with 
the state (blue) and site (red) return probability, we present the case 
of the 2-dimensional absorbing subspace span{^(l, 0, ^), (0, 1, 0)} 
(green). The upper and lower figure correspond to the Grover coin, 
Cg, and the Fourier coin, Cf^, respectively. 



Lattice 


Coin 


Approximate eigenvalues of II 


Square 




0.6593, 0.4069, 0.4069, 0.2878 


Square 


c% 


0.5517, 0.3882, 0.3882, 0.2880 


Hexagonal 




0.8017, 0.2411, 0.2411 


Hexagonal 


CI 


0.6365, 0.6365, 0.5462 



Here Cq is the three-dimensional coin 

1 /O y2\ 

This points to a general property whose explanation and implications deserve future 
research. 
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Appendix A. The average of the expected \^-return time 

Lemma A.l. If dim < oo and U has no eigenvectors in V-^, the eigenvalues of 
U = {I — P)U and A = {I — P)U{I — P) lie inside the unit disk. Hence, \\U'^\\f (ind 
WA^Wf converge to zero exponentially fast as n ^ oo. 

Proof. Let V G H \ {0} such that U%1) = Xtp, X e C. If Ui/j G V^, then Ui; = Ui; = 
Xip. Hence, if no eigenvector of U belongs to V^-*", we have Uip ^ V'^. But then 
lAlllV'lHr^lKllf/V^II = 11^11, so |A|<1. 

Consider now (j) El-iX {0} such that Acf) = X(f), A G C. Then, = {I — P)(j) satisfies 
Uip = {I — P)A(f) = Xip. U ip = 0, then A(j) = 0, i.e. A = 0. Otherwise, we are in the 
previous situation that leads to |A| < 1. 

The asymptotics of ||f/"||F and H^^Hf follow directly from the previous results and 
the Jordan form of U and A. □ 

A finite average of the expected ^-return time requires a finitely supported spectral 
measure fj,{dX) of V. We know that the mass points of /x(c/A) must be eigenvalues 
of U and the mass at A is /x({A}) = PE{{X})P, where E{{X}) is the orthogonal 
projector onto the eigenspace associated with A. 

Theorem A. 2. IfV has a finitely supported spectral measure fi{dX), then j^^^ t{'iI)) dip = 
K/ dim V with K a positive integer that can he computed equivalently as 

K = ^dim(EfcF) = 5];rank/x({Afc}), 

k k 

where {Xk} are the mass points of fi{dX) and Ek = E{{Xk}) are the orthogonal pro- 
jectors onto the corresponding eigenspaces of U . 

Proof. Under the hypothesis, V lies on a finite-dimensional [/-invariant subspace "Ho = 
®kEkV . Due to the unitarity of f/, the orthogonal subspace l-Li = T-Lq is also U- 
invariant and the unitary evolution splits as f/ = Uo(BUi with Ui = U \ "Hj. Therefore, 
to compute the average of the expected V^-return time, we can use the operator Uq 
instead of U since no state of V escapes from T-Lo during the evolution. In other words, 
= as operators on V. 

So, in what follows we can assume without loss that = (BkE^V. Then, the 
eigenspaces of T-L are E^V, which implies that U has no eigenvectors in V-'-. Indeed, 
any eigenvector will have the form tp = Ek(p, (p & V, and the orthogonality require- 
ment ip G V-^ forces \\ip\\^ = {(p\Ek(p) = {(pl'ip) = 0. This means that we can use the 
results of the previous lemma. 

We know that K = X]n>i '^II'^^IIf- Since a„ = Utj^'^P — JJ'^P, we have lla^HI^ = 
||fjn-ip||2^ — ||[/'^P|||,. These identities combined with lim„_>.oo ?^||[/'^|||' = from the 
lemma, lead to = X]n>o II^'^-^IIf' "which can also be expressed as 

K =\\P\\l + J2\\A"'B\\l, A = {I - P)U{I - P), B = {I-P)UP. 

n>0 

Using the equality AA^ + BB'^ = I — P, we can write 

= TiiA'^BB^A^)") = - A^ = I-P. 
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Thus, ^„>o 1 1 ^"-B III' becomes a telescoping series whose convergence follows from the 
lemma, giving finally 

K = \\P\\l + ||/ - P||| = dim 7^ = dim{EkV). 

k 

To finish the proof simply note that kerP fl E^V = V-^ fl E^V = {0}. This proves 
that dim{EkV) = dim{PEkV) = rank(PEfeP) = rank n{{Xk})- □ 

The subspace "Hq = ®kEkV appearing in the proof of this theorem is the minimal 
[/-invariant subspace containing V. The numerator K is its dimension. 

In the case of a finite-dimensional state space "H we can also express K = dim?/ — !/, 
where u is the number of linearly independent eigenvectors of U in V-^. To see this, 
note that we can consider the points running over all the eigenvalues oi U - i.e., 
over all the mass points of E{dX) - bearing in mind that fj,{{Xk}) = PEkP will vanish 
when Afc is not a mass point of ^{d\). Then, denoting by = EkH the eigenspace 
of U associated with A^, it is clear that 

z/ = ^dim(yfcny^). 

k 

A simple calculation shows that {EkV)-^^'" = Vk P^V-^, which means that = E^V ® 
(Vfc n V^-*") is an orthogonal decomposition of the eigenspace Vk- Hence, dimV^ = 
dimEkV + dim(Vfc fl V-^). If dim'H < oo the eigenvectors of the unitary operator U 
span thus, dim?^ = dim and 

K = ^dim(Pfc1/) = dim?/ - 

k 

The reason we can use Lemma lA.ll in the proof of the previous theorem - even if 
dim = oo - is that the lemma can be applied to the minimal [/-invariant subspace 
Hq containing V whenever dim Ho < c>o, i.e. if all the states of V have a finite expected 
^-return time. In that case, the unitary operator Uq = U \ T-Lq and the orthogonal 
projector Pq onto 1-Lq satisfy 

[/"Po = ((Po - P)U,T = and A^Po = ((Pq - P)[/o(Po - PT = ^o- 

Thus, we find an exponentially fast convergence ||[/"Po||f, ||^"-Po||f — ^ 0, which im- 
plies that ||[/"P||i? — exponentially fast as well. Nevertheless, Lemma [A. II does not 
ensure these asymptotics for states with a finite expected ^-return time that do not 
cover all of Sy- Even in this case, the following result shows that a weaker asymptotic 
result holds. 

Lemma A. 3. Any V -recurrent state with r(^/;) < oo satisfies 

lim n||[/>|p = 0. 

n—>oo 

Proof. For any G Sy, the identity 

r = t{4)) = ^ n{sn-i - Sn) 

n>l 



|a„^||2 = ||[/"-i^||2 _ ||[/'^^/;||2 yields 

y^Sk-nSn], s„ = ||[/"Vir 

fc=0 
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On the other hand, the l^-recurrence of '0 can be characterized in terms of the non- 
increasing sequence s„ of V^-survival probabihties in any of the following equivalent 
ways 



n>l n>l 



OO. 



Therefore, if is V^-recurrent and r < oo the series Xln>o must converge. Oth- 
erwise, nSn — > oo and thus, for some index Uq, 

n>l ^ ' n>no 

Hence, J2n>o = t' < oo and nSn r' — r. Finally, r' — r 7^ leads to s„ ~ 
(r' — r)/n, in contradiction with the convergence of ^„>o Sn- D 

Appendix B. Site Schur functions for one-dimensional coined walks 

The purpose of this appendix is to calculate the Schur function f^{z) of a site 
subspace = t), l^^? i)} for a one-dimensional coined walk with arbitrary coins. 
We will perform this calculation first for a semi-infinite lattice, extending it afterwards 
to a finite lattice and finally to a doubly-infinite lattice. 

On the half-line, fx{z) is the Schur function of the matrix measure 



^^■^''^ - U.(A) WA) |X,.,,(A)P ) 

where Xf:{z) are the orthogonal Laurent polynomials generated by the CMV matrix 
with Verblunsky parameters (70, 0, 72, 0, ... ) and ^{dX) is the corresponding orthog- 
onality measure. 

The above matrix measure can be rewritten in terms of the orthogonal polynomials 
(pk{z) with respect to fi{d\) and their reversed ones (pl{z) = z^Tp^{z~^) using the 
relations [m |2S1 EH EHl |34] 

which yield 

The first step is to compute the Caratheodory function Fx{z) of fj,^{dX), 

TP (A= ( ^2x(2;) G2x{,z) 

\G2.{z) F2x+iiz 
where Fk{z) is the Caratheodory function of \(pk{z)\'^fi{dX) and 

/X + z 
J— A~Vfc(A) <^fc+i(A) i^idX), 

^^(^)= / ^A-'=-Vl.(A)^I.+i(A)/i(rfA). 
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Let f(z) be a Schur function with Schur parameters (7a;)- Khrushchev's formula 
states that [T8] 



FJz) 



1 + zf 



k-l( ^ 



)fk{z) 



zf^~\z)Uzy 



We will use the notation F{z) = Fo{z), which is the Caratheodory function related 
to f{z). 

On the other hand, Gk{z) and Gk{z) can be expressed in terms of the second 
kind polynomials Qk{z), generated by the Verblunsky parameters (— 7^), due to the 
relations 



X + z 
X-z 
X + z 
X-z 



j-fk{X) -ipk{z) 



fi{dX) = Qk{z) 



X^ 



^IW-vliz) 



fi{dX) 



J =0,1, 
J = 1,2, 



,fc-l. 



k. 



These relations lead to the equalities 

Gk{z) = z-^Lpk{z) [i^k+i{.z)F{z) + Vtk+i{.z)] , 
Gk{z) = z-'-'^liz) [vU,{z)F{z) - ni^,{z)] . 

The Verblunsky parameters provide a forward and a backward recurrence relation 
involving the orthogonal polynomials ipk{z) and their reversed ones [T6l [26l 128], 

Pk<^k+i{z) = z^k{z) - ^k¥>*k{z), PkZ^k{z) = <^k+i{z) + 7fc<^fc+i(^), 
pk<^l+i{z) = --fkZ<^k{z) + <^*k{z), pk<^l{z) = lk<^k+i{z) + <^*k+i{z), 

with analogous formulas for Qk{z), replacing 7^ by — 7fc. From these recurrence rela- 
tions we get by induction an expression for the iterates and inverse iterates of f{z), 
namely 

¥>k+i{z) 



^ l ^l{z)F{z)-ni{z) 

' zMz)F{z) + nk{zy 

Therefore, Gk{z) and Gk{z) are related by 

Gkiz) fk+iiz] 



Gk(z) 



f 



k-ll 



A simple connection among the functions Fk{z), Fk+i, Gk{z) and Gk{z) results 
from rewriting 



Fkiz) 



X-z 



Using this expression together with the recurrence for the orthogonal polynomials 
yields 

Pk{Fk{z) - Fk+i{z)) = -fkGk{z)+^kGk{z). 
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Finally, the above identity combined with the relation between Gk{z) and Gk{z) leads 
to 



Bearing in mind Khrushchev's formula, it is obvious that, as in the case of Fk{z), we 
can express Gk{z) and Gk{z) in terms of iterates and inverse iterates of f{z). Indeed, 
such an expression can be simplified using the Schur algorithm to relate different 
iterates or different inverse iterates, giving 

Gkiz) '^^'f'-'^'^ 



Gk(z) 



{l-^kzf'-'{z)){l-zp{z)fk+,{z)y 

2pkzfk+i{z) 

il-^kzp-^iz)){l-zf'^{z)fk+,{z)y 



We have rewritten the coefficients of the matrix Caratheodory function Frc{z) in 
terms of iterates and inverse iterates of the scalar Schur function f{z). Now it only 
remains to perform a tedious calculation using again the Schur algorithm for these 
iterates to arrive at 



Uz) = z~\F,{z) - h){F,{z) + h 



,-1 



l2xf (z) P2xf2x+l{,z) 
p2xP''~\z) -l2xf2x+l{z] 



In the case of a finite lattice of sites, the only difference from the semi-infinite 
lattice is that the measure fi{d\) is supported on 2N points and thus there exists only 
a finite segment of 2N orthogonal polynomials and Laurent polynomials. The only 
effect of this finiteness in the previous arguments is that the Schur function f[z) has 
a finite number of Schur parameters (70, . . . ,'y2N-2, 1)? and similarly for its iterates. 
Hence, the site Schur function has the same form as in the half-line case. 

As for a doubly-infinite lattice, we can obtain the corresponding site Schur functions 
by a limiting process on semi-infinite lattices. Consider a coined walk on the line 
with coins given by a doubly-infinite sequence of parameters {■y2x)xez- The transition 
matrix in the basis states is then a doubly-infinite CMV matrix with Verblunsky 
parameters (. . . , 0, 7.2, 0, 70, 0, 72, 0, ... ). 

li U = SC is the corresponding unitary step, let us define for each site Xq G Z 
a new unitary = © U^^ where is the identity on the strict left subspace 
®x<xoVx, and U^^ = S^^C acts on the right subspace (Bx>xoVx with the same coin 
operator as U, but a minimally modified shift 5*^^ to impose a reflecting boundary 
condition from the right at site Xq: S^Jxq.I) = |xo,t)- other words, Ux^■, splits 
into the trivial evolution on the strict left subspace and a coined walk on a half-line 
- the sites x > Xq - whose unitary step acts exactly as U on the strict right subspace 
®x>xoVx, although not on the subspace V^q itself. Clearly, converges strongly to 
U as Xq ^ —00, i.e. \imxo^^oo{U — Uxq)^^ = for any state ip. 
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Additionally, the analysis of the half-line shows that, under the Uxq evolution, the 
matrix Schur function of a site x > xq is given by 

f _ /72x/Jo"H^) P2xf2x+l{z) 

jx^A^) \p,ji^^-\z) -%j2x+iiz: 

with fk{z) and f^^{z) given by the Schur parameters 

(7fc,7fc+i,7fc+2,---) and (-7^, -7;^-!. • • • > -72a;o' 1)' 

respectively. 

The strong convergence U^q — )■ U implies the weak convergence of the spectral 
measures ii^^^^dX) — fi^^dX) corresponding to any site subspace Vx, which also 
implies the uniform convergence f^^^xi^) ~^ fxi^) compact subsets of the unit disk 
for the related Schur functions. Hence, 



72xf (z) P2xf2x+l{z) 
p2xP''~\z) -l2xf2x+l{,z] 



with f {z) = lim^.Q^_oo fx^ {z). As a limit of Schur functions under the topology 
of uniform convergence on compact sets, p^~^{z) is a Schur function too and its 
Schur parameters are the limits of the Schur parameters of Taking the limit 

Xq -00 in the Schur parameters (0, -72^.2, 0, -73^.4, -73^^, 1) of fl^-^{z) 
shows that p^~^{z) is characterized by the infinite sequence of Schur parameters 



(0,-72a;-2,0,-72x-4, •••)• 
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